TRI-PP-99-027 



Singlet Interacting Neutrinos in the Extended Zee Model and 

Solar Neutrino Transformation 



G. C. McLaughlinQ and J. N. NgQ 
TRIUMF, 4004 Wesbrook Mall, Vancouver, BC, Canada V6T 2A3 



Abstract 

We study the impact of Standard Model singlet neutrinos on neutrino flavor 
transformation. We focus on an extension of the Zee model which includes singlet 
neutrinos, and find that the best limits on the interactions of the singlet neutrinos 
come from astrophysical phenomena. Singlet neutrino-electron scattering will im- 
pact both the mattter enhanced flavor transformation potential as well as detector 
cross sections. If electron neutrino - singlet neutrino oscillations are responsible 
for the solar neutrino anomaly, then the limit on the singlet neutrino interaction 
strength is of order of the weak interaction scale. Zee model modification of Ur — e 
scattering also impacts solar neutrino transformation, although this interaction is 
more tightly constrained. 
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Recent data on atmospheric, solar and accelerator neutrino experiments have in- 
dicated that neutrinos have a finite mass and they exhibit fiavor oscillation phenomena. 
Although the data has yet to be confirmed they point to the possibility that the three 
known active neutrinos are insufficient to account for all observations. The mass squared 
differences required to fft the data are: 10~^^ eV^ < (5m^ < 10~^ eV^ , 6m'^ ~ 10^^ eV^ and 
0.2 eV^ < 

"^"^LSND — 2eV^ where the subscripts s, and a refer to solar and atmospheric 
oscillations respectively. At least one light standard model singlet fermion is required to 
reconcile the data with the standard two by two neutrino mixing explanations. 

In most studies of a SM singlet neutrino one assumes it has no interactions other 
than gravitation at low energies and its role is to supply a mass for the active neutrinos. 
An example of a singlet neutrino is a right-handed neutrino of the Dirac type. Another 
example comes from SO(IO) GUT model where the heavy SM singlet is a Majorana 
neutrino. In this model, the heavy neutrino has a mass of 10^^"^'' GeV and gives the 
active neutrinos a small mass via the seesaw mechanism. Its large mass prevents it from 
having direct experimental consequences at current available accelerator energies. 

Recently there has been renewed interest in models of radiatively generated neutrino 
masses with new physics occurring at or not much above the weak scale. An example is 
the Zee model which generates neutrino mass with the same mechanism that produces 
lepton number violation 0. The crucial ingredient is a lepton ffavor changing SU(2) 
singlet scalar with U(l) hypercharge Y=2. It is straight forward to include a light SM 
neutrino in the model [Q, ^, so as to accommodate all of the anomalous neutrino data. 
We observe that now this light singlet neutrino now has an interaction with a strength 
determined by the mass of the Zee scalar and its Yukawa couplings to the SM leptons. 
A detailed phenomenological study of the current limits on these parameters are given in 
Ref i. 

This paper is concerned with the effects of SM singlet neutrinos that have interactions 
as weak as or weaker than the normal weak interactions. These interactions can arise 
from a spin 1 particle exchange as in extended gauge models or a spin boson exchange 
such as the charged scalar in the Zee model. The interactions may be weak due to small 
couplings and/ or a large mass of the mediating particle. One advantage of the Zee model 
is that it has relatively few parameters. Most are constrained by terrestrial experiments. 
This makes it a economical model for the study of singlet interacting neutrinos (SINs) on 
astrophysical phenomenon. Since there is now a large body of solar neutrino data with 
more to come, we look into the effects of SINs on some of the proposed solutions to the 
solar neutrino problem. For the purpose of this paper we can ignore the mixing of the 
Zee scalar with the requisite two Higgs doublets without loss of generality. 

SINs can affect the study of the neutrino fluxes from the sun in two ways. They can 
impact on matter enhanced flavor transformation and secondly they alter the neutrino 
detector cross sections since they now interact with electrons albeit very weakly. The 



1 



interactions of the singlet will come into play in any scenario which involves the Mikheyev- 
Smirnov-Wolfenstein (MSW) mechanism of neutrino transformation between an active 
neutrino and a singlet. We illustrate both effects with the small angle solution to the solar 
neutrino problem. If electron neutrinos transform to singlet neutrinos in the sun, the 
singlet neutrino coupling in the Zee model is most tightly constrained by the neutrino- 
electron scattering data at SuperKamiokande 0. We also show that a model which 
produces neutrino mass, such as the Zee model, may provide additional interactions for the 
active neutrinos which will in principle occur in active-active transformation scenarios and 
detector cross sections. As we will demonstrate however, these parameters are constrained 
by terrestrial experiments to have a much smaller effect on neutrino flavor transformation. 

We begin by looking at the general equation governing the transformation of electron 
neutrinos into another type of neutrino, where x = /i,r, or s, and s is the singlet neutrino, 
in a matter environment. This is given by 



ih 



d_ 

dr 







V9(r) 






. VA + Ke(r) 



A + VUr) 



where 



K(r) = ±2V2Gf 



V^{r))E-5m^ cos 2^^ 



In these equations 



A 



N:{r)-Nt{r)- 



5w? . 
^sm20„ 



(1) 

(2) 
(3) 

(4) 



= m\ — m\ is the vacuum mass-squared splitting, 9^ is the vacuum mixing angle, Gp 
is the Fermi constant, and N~{r), N^{r), and Nn{r) are the number density of electrons, 
positrons, and neutrons respectively in the medium. In the formulas given here, upper 
signs (in this case plus) correspond to the mixing of neutrinos while the lower signs (in this 
case minus), correspond to the mixing of antineutrinos. The potential ^Px{r) will have a 
standard model value and an additional term due to the extra interactions introduced by 
the Zee model. Other models which produce neutrino masses such as R-parity violating 
supersymmetric models may have similar interactions. Focusing on the Zee model, the 
charged scalar, h~ , is constrained by experiments at LEP to have a mass [Mh > 100 GeV). 
This in turn induces the following four fermion effective Lagrangian: 
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which governs low energy z/^e scattering. Note that this term has the opposite sign as 
the SM charged current interaction and is a prediction of the modeL This term has to be 
added to the standard model muon neutrino MSW potential as follows: 



V^{r) = T2V2 Gi 



N (r) 
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where 5 = ^/2\fi2\'^ /{SMlGp). The coupling constant I/12P is constrained by the mea- 
surements of the lifetime of the muon to be I/12P/ ( loooev )^ ^ .0015 (see [^). Therefore, 
S < .002. Similarly, z/g ^ Ur oscillations will be influenced by the Zee model via additional 
Ut- — e scattering contribution to the matter potential. The potential takes the form of 
Eq. I with the replacement of by K and 6 = ^/2\fi^\'^ / {%MlGp). The limit on [/igp 
is derived from LEP and SLC [0 measurements of the leptonic vector and axial vector 
couplings in Z decay. For a scale mass of = 800 GeV; the limit is 5 < 0.1. Clearly the 
effect can be much larger for Ue — Vt oscillation. This upper limit gives rise to a maximum 
change in the MSW potential of about 10%. 

The extension of the Zee model which includes a singlet neutrino produces neutrino- 
electron scattering terms which, after Fierz transformation, have the form 

Here, gi is the Zee model coupling between the singlet neutrino, the right handed electron 
and the scalar. This produces a singlet neutrino MSW potential of 

K(r) = GFf3{N,-{r) - iV,+ (r)) (8) 

where 

The best terrestrial limit on j3 comes from the leptonic right handed coupling to the 
Z [|]. The singlet neutrino does not couple directly to the Z boson, it only makes a 
contribution to the decay through a correction at one loop order. The limit on |(7ip/2M| 
is about 2 x lO'^^GeV"^, therefore /3 < 6. Other than this one loop effect we found no 
direct experimental bound since this singlet neutrino can only enter in weak interaction 
processes via leptonic mixing and no usable constraint is available. This in principle 
allows the singlet-electron interaction to be larger than the weak interaction. However, 
as we shall see below astrophysical considerations can provide tighter constraints. 

It is a characteristic of the Zee model that there are no interactions which convert 
electron neutrinos directly to muon, tau or singlet neutrinos by way of electron scattering 
mediated by the charged scalar. Therefore, in this model V^x = V^e ~ 0. 
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However, there are neutrino electron scattering terms in the Zee model which convert 
muon (tau) neutrinos to singlet neutrinos and vice versa. These are governed by the 
Lagrangian, 

1^ = - ^^2M^^ {^R^Kr)eReL - ^T^R(T'"'iy^{r)eRa^^eL^ + h.c. (10) 

In the forward scattering direction, the scalar term is proportional to neutrino mass and 
can therefore be neglected. The tensor term is proportional to electron spin and integrates 
to zero for unpolarized electrons. Therefore, for most situations, V^s{ts) = Vg^^sr) = 0. 

Returning to the case of Uf. — Vg transition, the combined active-singlet neutrino trans- 
formation potentials can be cast in the form: 
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where Nj^ is the total number density of nucleons. The electron fraction is defined as 

Ne-(r)- Ne+(r) , , 

A non-zero /5 will have the effect of increasing the potential if the electron fraction is 
greater than 1/3. In the sun for example, the electron fraction ranges from a value of 
about two thirds at the center to more than 0.85 in the outer layers. The main effect of 
the new potential is to change the position at which a given neutrino undergoes the MSW 
resonance. The position of the resonance is determined by the condition: 

Ve-Vs = ^m^co^1Q^. (13) 

Therefore, for given mixing parameters and a constant electron fraction, increasing /3 
causes the resonance position for a neutrino of energy E to shift to lower density. 

This scenario may be applied to several phenomenon, such as solar neutrinos and 
supernova neutrinos for the v^, Vg or v^. ^ z/^ situation. For the up-down asymmetry 
in the atmospheric neutrino problem, the oscillations between muon neutrinos and either 
tau neutrinos or singlet may be analyzed in a similar manner, although we note that in 
general there will be off-diagonal terms for z/^ ^ mixing. 

Taking again the example of the sun, we see from Eq. |T^ that all neutrinos will pass 
through the resonance condition at a position that is further from the center of the sun 
as /5 increases. For fixed brr? and sin^ 2^„, low energy neutrinos that in the case of /? = 
did not encounter resonances in the sun, will do so now if /? > 0. We illustrate this point 
in Fig. m where the survival probability for solar neutrinos is plotted for two values of the 
parameter This figure was produced by numerically integrating Eq (|l|) for the singlet 
neutrino. For /? = the solution reduces to the small angle sterile neutrino oscillation 
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solution to the solar neutrino problem, as in given in, for example |^ . It can be seen that 
increasing f3 will cause a decrease in the number of low energy electron neutrinos coming 
from the sun. 

In contrast, the nonzero (3 has little effect on the high energy neutrinos. The effective 
weak potential scale height 

d\n{Ve-Vs] 



dr 
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which determines the survival probability at the resonance position remains fairly constant 
with small changes in (3. In fact for a fixed and an exponential density profile and a 
given neutrino energy, it can be shown that the weak potential scale height remains 
constant regardless of the value of (3. Plots similar to Figure |l| can be drawn for z/g ^ u^ .,- 
oscillations, although nonzero 6 will have a smaller impact on the survival probability as 
the constraints on 6 are tighter. 

It is also important to take into account the effect of a non-zero f3 in neutrino detectors. 
In the Zee model, the neutrinos have only additional interactions with other leptons and 
not with the quarks due to the weak charge of the scalar. Hence, the radiochemical solar 
neutrino experiments such as SAGE and GALLAX will not be impacted by the 

additional interactions the singlet neutrino has. However, they will register the low energy 
neutrino flux which depends on (3 as explained before. On the other hand, an experiment 
which detects neutrino-electron scattering, such as SuperKamiokande however, will have 
some portion of its signal coming from singlet neutrino-electron scattering if singlet neu- 
trinos are present. The cross section for singlet neutrino electron scattering at first order 
is given by: 

where T is the electron recoil energy. In comparison, the largest contribution to the 
standard model z/g — e scattering cross section approximately given by: 
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The theoretical rate per electron recoil energy can be calculated by folding the survival 
probability for the electron neutrinos and the oscillation probability for sterile neutrinos 
with the cross sections and the fluxes of neutrinos. The detector rate may be estimated 



by folding the theoretical rate with an energy resolution function as in Eq 4 of Ref JTl . 

If Uf, ^ Vg transformation is the solution to the solar neutrino problem, then it can 
be readily seen that a constraint on (3 comes from the overall number of events in Su- 
per kamiokande. Various combinations of mixing parameters and values of j3 will produce 
different total count rate. If any significant mixing of electron neutrinos and singlet neu- 



trinos takes place in sun, and the flux predicted by the standard solar model is correct 



then j3 must be smaller than ~ 1/2. This limit is much stronger than the best limit from 
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accelerator experiments oi (3 < 6 and is derived from the extreme situation where all 
electron neutrinos above 5 MeV are converted to singlet neutrinos. 

Figure |^ plots the ratio of rates with to without matter enhanced flavor transforma- 
tion, using the shape of the neutrino spectrum from Several curves are plotted 



representing several mixing parameters and values of (3. It is seen that increasing the 
value of P causes a flattening of the recoil spectrum curve. It does not account for the 
upturn in event rate at high energy observed in the SuperKamiokande data 0. 

Correctly reproducing the overall rates in all of the solar neutrino experiments with 
a singlet neutrino - electron scattering interaction, requires an adjustment to the usual 
sterile neutrino MSW mixing parameters. For example if jS = 0.3, then must be 
increased by ~ 30% (see Eq. in order to avoid reducing the fluxes of pp neutrinos. 
The mixing angle must be adjusted to take into account both the effect of the change 
in on the survival probabilities of the '^B neutrinos and the effect of the nonzero 
singlet neutrino- electron scattering cross section in Kamiokande and Superkamiokande. 
For P = 0.3, sin^26't, must be increased by ~ 7%. The survival probability and expected 
electron recoil spectrum for these parameters is shown by the dot-dashed lines in Figures 
|] and 1^. 

We turn to the case of z/^ — z/g oscillations. The z/^ — e scattering cross section can 
increase by a maximum of a factor of 2, depending on the neutrino and electron energy, 
since the Zee model amplitudes and the standard model amplitudes add coherently. The 
scattering cross section in this case looks like 



da G%'m 



dT 2tx 



2 / T \ 
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For 5 = 0, this reduces to the standard model neutral current cross section. We illustrate 
the effect for matter enhanced flavor transformation with parameters = 5.4 x 10~^ eV^ 
and sin^ 29^ = 6.3 x 10^^. In order to reproduce the observations in solar neutrino 
experiments with nonzero 6, the change in the MSW potential will force a maximum 
increase in (5m^ of 10% above the 6 = solution. For this change in (5m^ the mixing 
parameter, sin^ 29^ retains approximately its original value in order to take into account 
the both change in and the increase in scattering at the detector. 

For comparison we consider the case of vacuum oscillations. There is no change in the 
survival probability of electron neutrinos due to the singlet interactions. However, the 
electron recoil spectrum will be effected. The singlet neutrino - electron scattering gives a 
signal similar to the neutral current scattering, although the size of the effect depends on 
the unknown magnitude of (3. Figure ^ shows the electron recoil spectrum for a vacuum 
solution of = 6.6 x 10~^^ and sin^26'^ = 0.9. The active neutrino solution is very 
similar to the singlet neutrino solution with (3 = 0.3. Therefore in the vacuum case, 
as in the small angle case, it may be difficult to differentiate between singlet neutrino 
oscillations and active neutrino oscillations just by detecting neutrino-electron scattering. 
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On the other hand, SNO [15| may be able distinguish the active-singlet solution from 
both the active-active oscillation and the active-sterile oscillation through the combination 
of the three reactions: (1) Ue + d p+p + e~, (2) u^ + d p + n + u^ and neutrino-electron 
scattering (3) + e~ ^ + e~ . In reaction 2, can be z/g, i/^, and u^-. In reaction 3, 
z/j. includes the three active neutrinos and also the singlet neutrino. The expected event 
rate per kilotonne year from the three reactions may be found in [|15| . The standard solar 
model predicts 6500 events from reaction 1 above a 5 MeV threshold, and about 2200 will 
be seen for vacuum i/e ^ i't oscillations. These numbers remain roughly unchanged in the 
presence of a singlet or sterile neutrino, although there will be some variation depending 
on the choice of mixing parameters. The number of events for reaction two is 710 from 
the standard solar model and remains unchanged in the case of active-active oscillations, 
although it should be reduced by about a factor of two in the case of active-sterile or 
active-singlet oscillations, depending on the mixing parameters. From reaction 3, about 
320 events are expected for active-active oscillations with about 78 coming from z^^(t)- The 
number of neutrino-singlet interaction events depends on the strength of the interaction 
(3. Therefore, if reaction 2 indicates sterile or singlet neutrinos, for large (3 reaction 3 in 
combination with 1 must be used to distinguish the two cases. 

In conclusion, we find that the extended Zee model with its new lepton number violating 
interactions gives rise to new neutrino-electron scattering mechanisms. The additional 
scattering occurs for muon, tau and singlet type neutrinos. It alters the small angle 
MSW Vf. — Vg and Vf. — solar solution by shifting the resonance position for neutrinos 
of a given energy. Furthermore for it increases the number of counts in water detectors 
and modifies the shape of the electron recoil spectrum produced from neutrino electron 
scattering. Taking these into account we have obtained a limit on singlet neutrino - 
electron interaction (5 which is an order of magnitude better that derived from current 
accelerator experiments. With non-zero Zee model interactions for the singlet and tau 
neutrinos, the apparent bvn} as measured by solar neutrino experiments can differ by as 
much as ~ 50% and ~ 10% respectively, from that which would be measured by reactor 
experiments. 

This work is partially supported by a grant from the Natural Science and Engineering 
Council of Canada. 
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Figure 1: Survival probability for electron neutrinos produced near the center of the 
sun is plotted against neutrino energy. The solid line is for electron neutrino-singlet 
neutrino mixing parameters of 5w? = 4 x 10~^eV^, sin^ 26v = 0.01 and a singlet neutrino 
interaction parameter, f5 = 0. This corresponds to the case of sterile neutrinos. The 
dashed line is for Srn^ = 4 x 10~®eV^, sin^ 29^ = 0.01, P = 0.3 and shows that increasing 
the value of (3 causes less electron neutrinos to survive. The dot-dashed lines represents 
Sm'^ = 5.2 X 10-6eV^ sin^ 2^^ = 0.011, p = 0.3. The solid and dot-dashed lines will 
produce very similar signals in current solar neutrino detectors. 



Figure 2: The ratio of the number of electrons produced by neutrino electron scattering 
to the predicted number of electrons using an undistorted spectrum in a detector such 
as SuperKamiokande is shown. The solid, dashed and dot-dashed lines correspond to the 
same parameters as in Figure |l]. Increasing j3 causes a flattening of the curve and overall 
increase in the number of counts. 



Figure 3: Shows the same ratio is as in Figure ^ for vacuum oscillations. The solid 
line corresponds to vacuum oscillations into active neutrinos, where the dashed line cor- 
responds vacuum oscillations into singlet neutrinos with interaction parameter f3 = 0.3. 
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